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The magnetized MINOS Near Detector, at a depth of 225 meters of water equivalent (mwe), is 
used to measure the atmospheric muon charge ratio. The ratio of observed positive to negative 
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atmospheric muon rates, using 301 days of data, is measured to be 1.266 ± O.OOl^tat.^o'oJ^syst.). 
This measurement is consistent with previous results from other shallow underground detectors, 
and is 0.108 ±0.019(stat. + syst.) lower than the measurement at the functionally identical MINOS 
Far Detector at a depth of 2070 mwe. This increase in charge ratio as a function of depth is 
consistent with an increase in the fraction of muons arising from kaon decay for increasing muon 
surface energies. 



PACS numbers: 13.85.Tp. 95.55.Vj, 95.85.Ry 



I. INTRODUCTION 

High energy cosmic-ray primaries interact with nuclei 
in the upper atmosphere and produce showers which 
contain pions and kaons. These secondary mesons can 
either interact in further collisions in the atmosphere, 
or decay to produce atmospheric muons. Since the 
majority of primary cosmic-rays are protons, there is 
an excess of positively charged mesons in the showers, 
and consequently, the atmospheric muon charge ratio 
/Na- is larger than unity. A precise measurement of 
the atmospheric muon charge ratio can be used to con- 
strain cosmic-ray shower models and, since atmospheric 
neutrinos are produced in conjunction with atmospheric 
muons, better determine atmospheric neutrino fluxes. 



decay probability. The larger value of €k implies that 
the kaon contribution to the muon flux, and therefore to 
the charge ratio, will be more significant at values of the 
vertical muon surface energy E^cosfl* exceeding e^. The 
charge ratio is expected to increase as a consequence 
of the increasing kaon contribution to the muon flux 
because in high energy interactions, single A +, s can 
be produced in associated production with strange 
baryons while single A"~'s cannot. The muon charge 
ratio is expected to decrease at even higher energies 
due to heavy flavor production 0, 0] and changes in 
composition of the cosmic-ray primaries @43- However, 
as these latter two processes are expected to only affect 
the charge ratio at energies greater than those studied 
here they are not considered further. 



The differential muon production spectrum in exten- 
sive air showers can be parameterized as [l[ : 

dN^ ^ 0.14E~ 2 - 7 / 1.0 0.054 
dE^ ~ cm 2 s sr GeV * I i I hl<^l + i i j~jjvjjjg 

(1) 

where is the muon surface energy and 9* is the zenith 
angle at the muon production point. Accounting for the 
curvature of the earth, this angle is geometrically related 
to the observed zenith angle 8, 

where R e is the radius of the Earth, and h is the muon 
production height for horizontal muons. Assuming the 
mean column depth for the initial cosmic ray primary 
interaction to be 85g/cm 2 gives /i=30km[2|. The two 
terms in Eq. ([1]) represent the contribution to muon 
production from pion and kaon decay respectively. The 
values 6^=115 GeV and 6^=850 GeV are the critical 
energies at the muon production height above which 
the pion and kaon interaction probability exceeds the 
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Atmospheric muon charge ratio measurements uti- 
lizing deep underground detectors are typically higher 
in value than those produced using detectors which are 
shallower since they sample muons with a larger value of 
' E^cosfl*. A previous measurement of the muon charge 
ratio utilizing the MINOS Far Detector, at a depth 
of 2070 meters of water equivalent (mweV, reported a 
value of 1.374 ± 0.004(stat.)iHio(syst.) for surface 
energies greater than 1 TeV. Another measurement Q 
of the charge ratio at TeV energies, conducted using the 
OPERA experiment located at a depth of 3800 mwe, is 
in good agreement with the MINOS Far Detector result. 
Atmospheric muon charge ratio measurements per- 
formed by the L3+C Q3, Bess-TEV [HI and CMS [H[ 
collaborations, for muons with surface energy E M above 
10 GeV and below 300 GeV, are consistent with the 2001 
world average of 1.268±[0.008 + 0.0002 • ^] [U]. 

The MINOS Near and Far Detectors Q are large 
underground magnetic calorimeters at depths of 225 mwe 
and 2070 mwe respectively. The detectors are designed to 
study neutrino interactions from the Fermilab Neutrinos 
at the Main Injector (NuMI) beam [l5j], but also trigger 
on atmospheric muons between beam spills. The depths 
of the MINOS detectors are ideal for making a measure- 
ment of the muon charge ratio at values of E^cos^* in 
the region dominated by the pion contribution and the 
transition region where the kaon contribution becomes 
significant. In this paper we present a measurement of 
the atmospheric muon charge ratio using the MINOS 
Near Detector. This result is then compared to the 
same measurement performed using the MINOS Far 
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Detector |[. 

This paper is organized as follows. Sections IU and Hill 
describe the MINOS Near Detector and the atmospheric 
muon Monte Carlo simulation respectively. Section IIVI 
describes the criteria for selecting atmospheric muon 
tracks with correct charge sign identification. Section IVl 
outlines the technique used to combine the data col- 
lected from two detector magnetic field configurations to 
determine the atmospheric muon charge ratio. Section IVl 
also elaborates on the evaluation of the systematic 
uncertainties. Finally, Sec. I VII presents the results of 
this analysis. 



II. THE MINOS NEAR DETECTOR 

The MINOS Near Detector [3 at Fcrmilab is a 
magnetized-steel and scintillator sampling calorimeter. 
It is located 94 m underground, with an approxi- 
mately flat overburden of 225 mwe. Its dimensions 
are 3.8mx4.8mxl6.6m. The detector contains 282 
vertical steel planes, each 2.54 cm thick. The scintillator 
layers are composed of either 64 or 96, 4.1cm wide and 
lem thick strips which vary in length from 2.5 to 4 m. 
The scintillator strips are rotated by 90° with respect 
to the previous layer to allow for three dimensional 
track reconstruction. The first 120 steel planes each 
have a scintillator layer attached to their surface. The 
scintillator layer of every fifth steel plane contains 
96 strips, and is said to be fully instrumented. The 
remaining four planes in each set of five are partially 
instrumented and contain 64 strips. Of the last 162 steel 
planes, only every fifth plane has an attached scintilla- 
tor layer. These scintillator layers are fully instrumented. 

Scintillation light is collected by wavelength shifting 
fibers embedded in the scintillator strips. Each strip 
is coupled to a single pixel on a 64-pixel multi-anode 
photo-multiplier tube |16| by a clear fiber. Each PMT 
pixel is digitized continuously at a frequency of 53 MHz. 
The detector response for a candidate atmospheric 
muon event is recorded when either four strips in five 
sequential planes, or when strips from any 20 planes, 
register a signal above the 1/3 photo-electron dynode 
threshold within 151ns. 

The detector's approximately toroidal magnetic 
field [ill varies in strength from 2.1 T near the magnetic 
coil hole to 0.2 T near the periphery of the steel planes. 
The magnetic field can be oriented to focus cither 
northerly going \i~ or These magnetic field orien- 
tations will be referred to as "forward" and "reverse", 
respectively. The detector is oriented horizontally 
and its long axis points 26.5548° west of north. The 
curvature induced by the magnetic field together with 
the three dimensional track reconstruction allows the 



determination of the charge sign of a muon traversing 
the detector. 

Detailed information regarding the MINOS Near De- 
tector, the electronics and the data acquisition systems 
can be found in [H, [13, ■ 



III. SIMULATED ATMOSPHERIC MUONS 

A sample of 42xl0 6 simulated atmospheric muon 
events are used to evaluate the analysis sensitivities 
and efficiencies, as well as to help assess systematic 
errors. Surface level muon events were generated by the 
HEMAS [3 atmospheric cascade simulation. Muons 
expected to intersect with the detector are propagated 
through a GEANT3 [2(| model of the overburden and 
the MINOS Near Detector. 



IV. EVENT SELECTION 

The selection criteria have been chosen to optimize the 
event selection efficiency observed in the data and the 
charge sign identification purity obtained from the Monte 
Carlo simulations. The data sample for this analysis 
consists of atmospheric muons with well-reconstructed 
energy and charge sign. Pre-selections are made to reject 
events that are not consistent with the passage of an 
atmospheric muon through the detector. Further selec- 
tions are applied to ensure track reconstruction quality 
and good charge sign determination. A summary of the 
selection criteria, as well as the selection efficiencies and 
charge identification purities, can be found in Table [H 

This analysis uses data collected between 2006 and 
2009. During the data taking period, the direction 
of the magnetic field was periodically changed be- 
tween the two orientations referred to as forward and 
reverse magnetic field. The final reverse magnetic 
field exposure was 150.5 days. Equivalent exposures of 
forward magnetic field data were collected immediately 
following the reverse magnetic field data periods to 
reduce systematic uncertainties. A total of 7.16x10 s 
triggers were collected over a combined 301 day exposure. 



A. Pre-Selections and Track Quality 

For an event to be included in this analysis it must 
possess a single downward-going atmospheric muon 
track, have been collected during a period of good 
detector run conditions, and have a reconstructed initial 
interaction point (RIIP) within 50 cm of the detector 
edge. If the reconstructed initial interaction point is 
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further than 3 cm outside the detector volume the track 
is rejected. The curvature of the track is known to 
be poorly determined, and the track is rejected, if any 
scintillator strips hit are further than 3 cm from the 
reconstructed position of the track, or if the track does 
not pass through a region of the detector where there is 
scintillator on each layer of steel. 

Track reconstruction errors may occur when the 
event contains a large amount of activity which is 
not related to the track. These extra hits degrade 
the charge sign determination and can be generated 
by muon bremsstrahlung, natural radioactivity, or by 
electrical and optical cross-talk between the channels on 
the multi-anode PMT [l6|. Events are rejected if more 
than 40% of the strips hit are not track-related. Muon 
tracks determined to be poorly reconstructed by internal 
consistency checks of the reconstruction algorithm are 
also removed from the data sample. 



B. Charge Sign Quality Selection 

Two selection variables are used to further increase 
the degree of confidence in the assigned curvature and 
charge sign of the track. The Kalman filter [2lJ used 
in the track curvature fitting provides an uncertainty 
a(q/p) on the measured value of q/p, where q is the 
charge and p the momentum of the track. The first 
charge sign quality selection is based on the value of 
(q/p)/a(q/p) determined by the track fitter. The second 
selection variable BdL is defined to be equivalent to 

N 

I Bi x dLi\ where N is the total number of planes 

i=l 

in the muon track, and Bi and dLi are the magnetic 
field and muon path length vectors, respectively, at 
plane i. A selection based on BdL is used to ensure 
that the magnitude of the bending due to curvature 
is larger than the apparent bending due to multiple 
scattering. Figure [T] shows the muon charge ratio as a 
function of (q/p)/cr(q/p) and Fig. [5] the charge ratio as 
a function of BdL, for data collected in both magnetic 
field orientations as well as the combined data set. The 
charge ratio for data collected during a single magnetic 
field orientation is defined as the ratio of positive 
to negative muons collected only in that orientation. 
The observed variation in the charge ratio above the 
selection thresholds in a single field orientation, as 
well as the difference in the charge ratio between the 
two different field orientations, stems from acceptance 
effects due to the magnetic field, detector asymmetry 
and detector alignment errors. The technique used to 
combine data taken in the two field orientations and 
remove these biases is discussed in the next section. For 
this analysis we required that BdL > 3.0 T • m and 
(q/p)/a{q/p) > 3.0. Below these values the charge ratio 
tends towards unity, indicating a degradation of the 



charge sign determination. Events which have passed 
all the selections described in this section are used in 
the calculation of the atmospheric muon charge ratio 
described in the next section. 
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FIG. 1: Charge ratio as a function of (q/p)/u(q/p) after all 
selections and requiring that BdL > 3.0 T ■ m. The vertical 
line is the (q/p)/a(q/p) threshold value used in the charge 
sign quality selection. 



V. CHARGE RATIO DETERMINATION 

As demonstrated in the previous section acceptance ef- 
fects in the Near Detector introduce a bias in the charge 
ratio when it is calculated using only data from a sin- 
gle magnetic field orientation. Figure [3] shows the charge 
ratio as a function of azimuthal angle, a variable sensi- 
tive to these biases, for the forward and reverse magnetic 
field data sets. Canceling these biases is done in the same 
manner as described in 0, If e\ is the efficiency for 
the selection of p + and £2 is the selection efficiency of 
\x~ in the forward field direction (FF) then the selection 
efficiencies for p + and pT in the reverse field direction 
(RF) are £2 and £\ respectively. Two independent equa- 
tions for the charge ratio, r a and n,, can be constructed 
in which the acceptance effects cancel. These ratios, cor- 
rected for live time, are 

r« = W F /t FF )/(N^ F /t RF ), (3) 

and 

n = (K + F /tRF)/(Nf F /t FF ), (4) 
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Selection 


Data 
Selection 
Efficiency 


Monte- Carlo 
Charge ID 
Purity 


Number of Triggers 


7.f6xf0 8 (100%) 




Pre-Selections 




Single Track Events 
Detector Quality 


3.18 x10 s (44.46%) 
3.16 x10 s (44.19%) 


70.7% 
70.7% 


Track Quality Selections 




RIIP and Curvature Selections 
Track-related Activity 
Good Reconstruction 


1.01x10 s (14.19%) 
6.13xl0 7 (8.56%) 
5.31 xlO 7 (7.42%) 


72.9% 
77.5% 
78.8% 


Charge Sign Quality Selections 




(q/p)/a(q/p) > 3.0 
BdL >3.0 T-m 


l.llxlO 7 (1.55%) 
3.23xl0 6 (0.45%) 


97.0% 
99.5% 



TABLE I: Summary of the event selection. Each row shows the total number of events remaining after all the applied cuts; 
in parenthesis, the percentage of events remaining; and lastly, the percentage of the remaining Monte Carlo events which have 
their charge sign determined correctly. 
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FIG. 2: Charge ratio as a function of BdL, after all selections 
and requiring that (q/p) /a(q/p) > 3.0. The vertical line is the 
BdL threshold value used in the charge sign quality selection. 



where N 11 (N^ ) is the number of positive (negative) 
muons observed during an exposure time t in a particular 
magnetic field orientation. The geometric mean of r a and 

n 



N„- 



s/r a n 



\ 



(5) 



provides a measurement of the charge ratio that is free 
of biases due to geometric acceptance effects, alignment 
errors and the different magnetic field live times. 

Figure [3] illustrates that the significant bias in the 
charge ratio, that is apparent in a single field orientation 
data set, is strongly suppressed in the geometric mean. 

After the selection criteria in Table U have been 
applied we obtain a final data set of 3,234,066 events. 
In the forward field sample we select 893,854 fi+ and 
721,428 \x~ . In the reverse field sample we select 912,944 
/i + and 705,840 n~ . The resulting charge ratio obtained 
by applying Eq. © is N „+/JV /i -=1.266±0.001(stat.). 



A. Systematic Uncertainties 

The event selection criteria are chosen to remove 
events in which the charge sign of the muon track has 
been incorrectly assigned. Systematic uncertainties 
associated with event selection are determined by 
establishing the sensitivity of the measured charge ratio 
to variations in the selections above their thresholds. 
Two selections, track- related activity and (q/p)/<r(q/p), 
display variations in the observed charge ratio above 
their thresholds. These data were divided into three 
equal statistics samples, corresponding to increasing 
confidence in the muon charge sign, and their charge 
ratios calculated. The maximum deviation from the 
nominal charge ratio, 0.012 and 0.006 for the track 
related activity and (q/p)/a(q/p) selections respectively, 
is taken as the systematic uncertainty associated with 
that particular selection. 
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FIG. 3: The observed charge ratio as a function of azimuthal angle. The charge ratio varies as a function of azimuthal angle 
in the forward and reversed data sets due to acceptance effects and alignment errors. When data from the two fields are 
combined using Eq. ([5| a flat distribution is obtained, indicating that the residual uncertainty due to these biases is small. The 
uncertainties shown are statistical only. 



Another systematic uncertainty relates to the re- 
maining events that have a misidentified charge sign 
after the selections. Since the atmospheric muon charge 
ratio is a value greater than unity more positive than 
negative muons will have their charge sign misidentified. 
Thus charge sign misidentification can only decrease the 
measured charge ratio resulting in a positive one-sided 
systematic uncertainty. Monte Carlo studies suggest 
that 0.5% of the events in the final data sample have 
an incorrect charge sign determination. However, one 
cannot discount the possibility that the misidentification 
rate is different in the data than in the Monte Carlo 
simulation. For this analysis it is assumed that the 
true charge sign misidentification rate in data differs 
from that in the nominal Monte Carlo by an energy 
independent factor a. The true charge sign misiden- 
tification rate in data can be estimated by exploiting 
the fact that it influences the shape of the charge ratio 
versus (q/p)/a(q/p) curve shown in Fig. [T] A toy Monte 
Carlo was written to produce similar curves assuming a 
charge ratio R a , and a charge sign misidentification rate 
equivalent to that in the nominal Monte Carlo scaled by 
a factor a. The difference between the measured charge 
ratio and R a for the values of a and R a that provide the 
best agreement between the data and the Monte Carlo 
prediction is taken as a systematic uncertainty on the 
measured charge ratio. Figure 2] plots the charge ratio 
as a function of (q/p)/a(q/p), along with the nominal 
and best-fit Monte Carlo. The best fit Monte Carlo is 
obtained using a charge ratio R a of 1.272. A similar 
study was performed on the curvature of the charge 
ratio versus BdL data in Fig. [2j in that study the 



best fit Monte Carlo is obtained with an actual charge 
ratio of 1.266. The maximum deviation of R a from the 
nominal charge ratio is +0.006, which is taken as the 
one-sided systematic uncertainty associated with charge 
sign misidentification uncertainties. 

Two alternative charge confidence selection criteria, 
similar in motivation to the BdL, have been investigated. 
The first alternative required that the muon traverse 
at least 27 planes at positions covered by the partially 
instrumented scintillator planes, within 1.75 m of the 
magnetic coil hole. The second alternative required that 
at least 37 planes be traversed where the \B x dL\ for 
each of those planes is greater than 0.03 T-m. These 
particular selection criteria have been chosen to optimize 
data selection efficiency while maintaining a charge 
misidentification rate in the Monte Carlo similar to that 
of the default BdL selection. The maximum deviation 
from the nominal charge ratio was observed to be 
0.003. This value is taken as the systematic uncertainty 
associated with the BdL selection. 

An imperfect reversal of the detector's magnetic 
field would introduce a geometric bias and lead to a 
systematic uncertainty in the charge ratio stemming 
from imperfect acceptance asymmetry cancellation. 
This error can be determined by evaluating the charge 
ratios calculated using Eq. §5§ and Eq. since the 
quantities r a and r& will diverge as the magnitude of 
this bias increases. It is found that r a and r& agree 
to within their statistical uncertainties, indicating that 
the magnetic field bias is negligible compared to the 
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FIG. 4: Charge ratio as a function of (q/p) /a(q/p), compared 
with the nominal and best-fit charge misidentification Monte 
Carlo simulations. 



Uncertainty Classification 


A— ei 


Track-Related Activity 

(<?/p)Mg/p) 

BdL Selection 

Charge Misidentification 


±0.012 
±0.006 
±0.003 
±0.006 


Total Systematic Uncertainty 


+0.015 
-0.014 



TABLE II: Summary of the systematic uncertainties in this 
analysis. The total systematic uncertainty is the quadratic 
sum of the individual uncertainties The magnitude of these 
uncertainties are independent of E,j and E^cos^*. 

statistical uncertainty. 

Tabic [XT] lists the systematic errors that are dis- 
cussed above, and their assigned values. Including 
systematic uncertainties the atmospheric muon charge 
ratio measured at the MINOS Near Detector is 
1.266±0.001(stat.)l~°;°^(syst.). 



B. The Muon Charge Ratio Underground 

The maximum momentum for which the charge sign 
of a muon track can accurately be determined with 
the MINOS Near Detector is limited by the scintil- 
lator granularity and the strength of the magnetic 
field. The reconstructed momenta p^^et of all muon 
tracks which pass the selections outlined in Table U 



is shown in Fig. O The maximum momentum of 
muons used in this analysis is approximately 300 GeV; 
the mean momentum is 15 GeV. The observed charge 
ratio as a function of reconstructed track momentum 
underground is shown in Fig. [6] The charge ratio is 
observed to be flat as a function of track momentum 
to within the statistical uncertainty of the measurements. 




FIG. 5: The reconstructed momentum of the muon tracks 
which survive all the selection criteria outlined in Table [T] 
The mean momentum is 15 GeV. 



VI. THE ATMOSPHERIC MUON CHARGE 
RATIO AT THE SURFACE 

The muon energies measured at the Near Detector 
depth must now be converted to energies at the surface 
of the Earth by accounting for energy lost by the muons 
in the overburden above the detector. For muons, the 
energy loss in matter can be parameterized by 

3 

--^ = a(E ll ) + J2b n (EJE li , (6) 

n— 1 

where X is the slant depth; a is the ionization term 
and the b n account for the radiative energy loss from 
bremsstrahlung, pair production and photo-nuclear 
processes. These parameters have a mild energy de- 
pendence, whose values for standard rock can be found 
in [23| . Using the technique outlined in j24|, the energy 
lost by muons traversing the Near Detector overburden 
has been calculated as a function of reconstructed track 
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FIG. 6: Charge ratio as a function of reconstructed under- 
ground track momentum. The y-axis uncertainties are statis- 
tical. The x-axis uncertainties are the width of the momentum 
bins used; the data are plotted at the median momentum val- 
ues. The dotted horizontal line is the best fit charge ratio of 
1.266. 



energy and zenith angle. These energies have been cal- 
culated assuming a flat vertical overburden of 224.6 mwe, 
which is comprised of two distinct geological layers. 
Above the cavern hall lies 72.1m of Dolomite/Shale 
bedrock with a density of 2.41 g/cm 3 , followed by 22.2 m 
of Glacial till with a density of 2.29 g/cm 3 . It has been 
suggested in that at the MINOS Near Detector 
depth the slightly higher rate of energy loss of /Lt + 
over fi~ |25M30| could lead to a surface charge ratio 
that is slightly higher than that observed underground. 
However, as the magnitude of this effect is negligible 
when compared to the systematic uncertainties of this 
measurement we have assumed the same energy loss 
function for both charges. 

The surface energy resolution obtained with this 
extrapolation method is dependent on the accuracy of 
the overburden model and the muon pointing accuracy. 
Topographical measurements [3lj show that the surface 
altitude varies by no more than 13 m within 3 km of 
the detector. This contributes a 14% uncertainty on 
the surface energy estimation at all zenith angles. The 
zenith pointing resolution of the Near Detector has been 
determined to be 1.1±0.2° by measuring the zenith 
angle separation between the two muons in multi-muon 
cosmic-triggered events for which both muons pass the 
charge ratio selection criteria. The corresponding error 
on the surface energy is negligible for vertical muons 



increasing to 5% at 70°, 14% at 81°, 25% at 85° and 50% 
at 87°. 

The extrapolated muon surface energy distribution, 
E M , is plotted in Fig. [TJa). The muons which populate 
the high energy tail have a large zenith angle and 
thus pass through the largest amount of matter before 
reaching the detector. To first approximation the 
minimum energy a surface muon needs to reach the 
Near Detector is 52GeV/cos6>. The E^cosfl* distribution 
is much narrower than the E M distribution as the 
cos6** and cos# terms effectively cancel. Due to this 
cancellation the muon events which occupy the high 
energy tail in Fig. [TJb) are those which possess the 
highest reconstructed momentum at the detector. 

Table IIIII presents the measured atmospheric muon 
charge ratio, and E M cos6>* in equal bins of cos#*. The 
median surface energy increases as cos(9* decreases due 
to the increasingly large overburden. E^cosfl* increases 
more slowly than E M as the increase is due primarily to 
the detector's larger analyzable momentum range as a 
function of zenith angle [32[. The charge ratio is observed 
to be independent of cos#*. E p and E M cos6** to within 
the uncertainties of our measurements. The charge ratio 
measurement at E u =964GcV is consistent with the 



MUTRON spectrograph measurement [22j, but is lower 
than other TeV energy scale measurements [1, [^, |33| . 
This difference is due to the range of zenith angles sam- 
pled by each of the detectors and the cos#* dependency 
of Eq. ©. 

The "ttK" @, S] model is derived from the differen- 
tial muon production spectrum parameterization given 
in Eq. ([T]), and predicts that the muon charge ratio is 
only dependent on E^cos^*. The 7rK model is a qualita- 
tive model describing the relative contribution of pions 
and kaons to the atmospheric muon charge ratio. Fol- 
lowing the prescription in Q and defining / w and fx as 
the fraction of all decaying pions and kaons which decay 
with a detected the atmospheric muon charge ratio 
N^+fN^- is given in Eq. Q: 



J L 

\ 1 + 



cose*/e„ 



0.054X/K 
1 + XXEp cosO*/e K 



\ 1 + l.l£„ci 



0.054x(1-/k) 
1 + LIE,, cos0*/e K 



(7) 



The charge ratio as a function of E^cosfl* is plotted 
in Fig. [8] along with the results from L3+C [10l|[34l. 
Bess-TEV [U, UTAH [H, MUTRON Q, OPERA [gj , 
CMS H2 and the MINOS Far Detector 0, ||. In 
instances where the data were not already listed BjS ci 
function of E^cos^* we convolved the published cos# 
acceptances with the given E^ data, with assistance 
from the original authors in the case of L3+C (HJ . 
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cos6>* 


Surface Energy 


E^cosfl* 


Char 


ge Ratio 






[GeV] 


[GeV] 




Interval 


Median 


Interval 


Median 


Interval 


Median 


r 


±(stat.) 


>0.9 


0.914 


55-67 


61.5 


53.5-61.0 


56.5 


1.257 


0.069 


0.8-0.9 


0.824 


60-77 


70.5 


54.2-62.9 


58.0 


1.270 


0.013 


0.7-0.8 


0.735 


69-91 


80.5 


54.6-65.0 


59.2 


1.277 


0.005 


0.6-0.7 


0.642 


79-112 


94.5 


55.2-69.0 


60.8 


1.270 


0.003 


0.5-0.6 


0.549 


95-142 


115 


56.2-74.0 


62.7 


1.269 


0.003 


0.4-0.5 


0.452 


115-191 


145 


57.4-80.9 


64.8 


1.259 


0.003 


0.3-0.4 


0.356 


148-280 


195 


59.4-91.0 


68.0 


1.263 


0.003 


0.2-0.3 


0.262 


212-484 


290 


63-109 


75 


1.263 


0.005 


0.15-0.2 


0.181 


382-864 


513 


76-140 


92 


1.276 


0.013 


<0.15 


0.139 


697-3450 


964 


104-390 


134 


1.249 


0.032 



TABLE III: The charge ratio, E M and E^cosfl*, in equal bins of cos#*. The range of energies observed in each cos#* bin is 
comparable to the energy resolution expected from the surface extrapolation. The charge ratio is observed to be independent 
of cos#*, and E^cosO* to within the uncertainties of our measurements. 



A x 2 P er degree of freedom fit test to the irK model 
was performed over (//£•,/„■) space using only the charge 
ratio measurements from the MINOS Near and Far De- 
tectors. The x 2 minimum was found at 7^=0.55 and 
fx=0.70. These results are consistent with earlier fits 
by The best fit curve to the data is plotted in 

Fig. [5] and indicates that the kaon contribution to the 
charge ratio becomes significant at E^cosfl* greater than 
a few hundred GeV. Using two functionally identical de- 
tectors, at two different depths, we have demonstrated 
that the increase in charge ratio observed at the deeper 
Far Detector is consistent with an increase in the frac- 
tion of observed muons arising from kaon decays in the 
extensive air shower. 



extrapolate the muon surface energy. No statistically 
significant change was observed in the charge ratio as a 
function of either the surface energy E^, from 62 GeV to 
960 GeV, or in the vertical muon surface energy E M cos#* 
from 57 GeV to 134 GeV. This work presents the first 
demonstration of an increase in the observed muon 
charge ratio between functionally equivalent shallow and 
deep underground detectors. The most likely source of 
this increase is the greater probability that a muon in the 
deep detector data sample results from kaon production 
in extensive air showers. A fit has been performed to a 
simple parametric model of muon production from pion 
and kaon parents, and the results of the fit support this 
interpretation of the combined MINOS Near and Ear 
Detector charge ratio measurements. 



VII. SUMMARY 



A charge ratio measurement has been performed on 
301 days of atmospheric muon data collected using the 
MINOS Near Detector. The atmospheric muon charge 
ratio measured at 225 mwe underground is 

— ^ = 1.266 ± 0.Q01(stat.)±%%\l(syst.). (8) 

The reconstructed underground energy and zenith 
angle for each muon used in the analysis were used to 
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FIG. 8: The atmospheric muon charge ratio as a function of E,jCos#*. The y-axis uncertainties are the statistical and systematic 
uncertainties added in quadrature. The MINOS Near Detector data, re-binned in five equal cos#* intervals, are plotted for 
each bin at the median E M cos#* values; the x-axis uncertainties are the bin widths. The dashed line is the best fit curve to the 
7rK model using only MINOS Near and Far Detector data. 



